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Abstract 

^ . We present exact analytical bright and dark (black and grey) solitary wave solu- 

^ \ tions of a nonlinear Schrodinger-type equation describing the propagation of spatial 

beams in media exhibiting a saturable nonlinearity (such as centrosymmetric pho- 
torefractive materials). A qualitative study of the stationary equation is carried out 
together with a discussion of the stability of the solutions. 
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Q , 1 Introduction 

K> . 

'jjj . The Nonlinear Schrodinger Equation (NLSE) in its many versions is one of the 

5^ \ most important models in mathematical physics, with applications to numer- 

ous fields [Tp] such as, for example, in semiconductor electronics [31H] . pho- 
tonics [5p] . plasma physics [7], fundamentals of quantum mechanics [S], dy- 
namics of accelerators [5] , mean- field theory of Bose- Einstein condensates [TU] 
or biomolecule dynamics [11], to name only a few of them. In some of these 
fields and many others, the NLSE appears as an asymptotic limit for a slowly 
varying dispersive wave envelope evolving in a nonlinear medium [12]. The 
study of these equations has served as a catalyzer for the development of new 
ideas or even mathematical concepts such as solitons [I3] or singularities in 
partial differential equations [T^IITS] . 

Among the many different variants of the nonlinear Schrodinger equation, one 
which has attracted a lot of interest is the so-called local saturable nonlinear 
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Schrodinger equation, which we write in the form 



.dA . Id^A . _„ .,2, 



'aF + ia^ + ^'^d-^l)'^ 



where A is a complex envelope function, and crV"(|Ap) represents the nonlin- 
earity. Equation ([1]) describes the propagation of a laser beam in a saturable 
nonlinear medium in (1+1) dimensions. A{x, z) measures the (complex) beam 
amplitude with x and z being the transverse and longitudinal coordinates, 
respectively. Here, I^d^p) includes the light-induced refractive index change 
and (T G M is proportional to the strength of the nonlinearity. For local sat- 
urable nonlinearities, V^(|Ap) oc (1 + lAp)""", with a > being the saturation 
index. Equation ([1]) has shown a number of applications in plasma physics [16] 
and nonlinear optics pnil8fl9] , specifically in photorefractive media. Photore- 
fractive materials manifest a wealth of nonlinear phenomena which include the 
propagation of solitons [20|21ll22ll23|24j . surface waves [26] and slow light [25], 
pattern formation [27], charge singularities [28], and critical enhancement |29j . 
The case with a = 1 (which corresponds to non-centrosymmetric photorefrac- 
tive media) has been throughly studied, and solitary wave solutions in the form 
of bright, black and grey solitons have been found numerically [201[^T][?2 ,30j. 
It was also shown theoretically (by means of numerical calculations) [21] , and 
experimentally [32][33] . that the nonlinearity with a = 2 also appears in the, so- 
called, centrosymmetric photorefractive media and exhibits bright and black 
spatial solitons. Motivated by these studies, here we go beyond the previous 
numerical studies and obtain three classes of exact solutions to Eq. ([T]) in 
closed form: bright, black and grey solitons (the latter class being addressed 
here for the first time, to the best of our knowledge). 

Although many methods have been used to construct exact solutions to various 
types of nonlinear Schrodinger equations such as the cubic one in different sce- 
narios [3^35|36|l37] . the cubic-quintic [38ii39j and others [IQ], it is remarkable 
that no exact solutions are available for the saturable nonlinear Schrodinger 
equation (^ ioi a = 2 despite its practical applications. 

The paper is organized as follows. In section 2, we first impose the necessary 
conditions to derive solitary waves of a general version of the nonlinear paraxial 
wave equation in which the nonlinearity depends locally on the light intensity 
(the latter being proportional to \A\'^). Section 3 presents a qualitative study 
of the stationary saturable nonlinear Schrodinger equation. Sections 4, and 5, 
contain the main results about the study of analytical solutions to Eq. ([1]): we 
obtain exact bright, and dark solitary waves, respectively, together with their 
existence curves, that constrain the amplitude and the width of the waves. 
Stability of bright and dark solutions is also discussed. 



2 Solitary Waves of the Nonlinear Paraxial Wave Equation 

For completeness, our first purpose is to deal with a general local nonlinearity 
cry(|y4p) to derive the characteristic wave equation for solitary waves from 
Eq. ([T]). It is more convenient to express the solutions of Eq. ([T]) in terms of 
the real amplitude u{x, z) and phase 0(x, z) as 

A{x,z)=u{x,z)e''^^''^'\ (2) 

Substitution of Eq. ([2]) into ([1]) yields the system of equations 
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It is clear that solving the coupled system of Eqs. (I3aj) and (l3bl) is equivalent 
to solving Eq. ([T]). 

We now impose the fundamental property of solitary waves: their shape in- 
variance with respect to Galilean boosts 



u{x, z) = u{x — vz) 



(4) 



where v can be regarded as the dimensionless transverse spatial velocity (or the 
steering angle between the propagation direction and the z— axis). In particu- 
lar, when f = 0, the profile of solitary waves remains unchanged with respect 
to translations along the z— axis. Let ^ = x — vz and C = -2- Then, the system 
of Eqs. fl3a|) and fl3bl) reads as 
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Equation fl5aj) can be readily integrated with respect to ^, and yields 

r« dr 



0(e,C) = C7(C) + < + ^(C)/ 



u^{t) ' 



(6) 



where C{() and -D(C) are two integration functions that depend solely on (. 
Also, since 0(^, () satisfies Eq. (I5b1l . inserting expression ([61) into Eq. (l5bl) one 



finds 



2k 



d_ 



C{0 + D{0f 



i dr 



u^{t) 



u^ u at, 



that must be fulfilled for all ^ and (. Now, if u{C,) is a non-constant function, 
this implies that the only admissible solutions for C{() and -D(C) are 



D{0 



^roC + 

■ ^0 ) 



(8a) 
(8b) 



where Fq is a constant (propagation constant), 00 is a fixed reference phase, 
and Iq is a parameter. Under these conditions, Eq. ([7]) reduces to the following 
characteristic equation for solitary waves 



d?u 

de 



^ - 2rM - 4 + 2aViu^)u = 



(9) 



with r = To — Y- -'^f ^0 7^ 0, in order to obtain nonsingular solutions to Eq. (jH]), 
it is necessary that m 7^ for —00 < ,^ < 00. Multiplying Eq. Qj by du/dC, 
and upon integration, it follows that 




Tu^ + ^ + aU{u)=So, 



(10) 



where U{u) = 2 /" V{T'^)TdT and Sq is a constant (it plays the role of an effec- 
tive total mechanical energy). It is expected that if the nonlinear contribution 
is absent, i.e. if aU{u) = 0, then, the solution to Eq. fITOl) . given by 



u\0 



,±V2r(?-Co) _ /i c e=Fv^(?-5o) 



+ /^e^ 



cannot represent a localized wave for any value of F, /q, and Sq. The presence 
of aU{u) is thus essential for the existence of solitary waves. Henceforth, we 
assume a 7^ 0. Moreover, for solitary waves to exist it is clearly required that 
£0 - aU{u) + Tu"^ - ll/{2u^) > 0. 

Since spatial solitary waves are well-localized excitations in space, we look 
for solutions such that their derivatives of all order vanish at infinity. Three 
families of solutions are distinguished according to additional boundary con- 
ditions at infinity and at the origin ^ = ^o- Let u = u^o for |^| -^ cxo, and 
u = Uq, ^ = u'q both at ^ = ^0- Solitary waves whose asymptotics tend to 
zero (Moo = 0), together with u'q = 0, are the well-known bright solitary waves. 
Those that fulfill Uq = (with u^o 7^ and u'q 7^ 0) are the black solitary 
waves, whereas those that do not vanish in — 00 < .^ < 00 and satisfy Uq = 



are grey solitary waves (both black and grey are denoted as dark solitons). In 
contrast with bright and black solutions, grey solitary waves possess nonzero 
values of Iq. Using the above prescribed boundary conditions, formulae for the 
constants F, Zq, and £q can be found for each family of solutions. 

Bright solitary waves correspond to < m^ < Mq, 
^ a [U{u,) - f/(0)] 



ul 



/o = 0, ^o = ^f/(0), (11) 



and are the implicit solutions obtained after integration of 



2 



2 



L(^A ^^(o)-t;(,.) + l^"'°'-f°)"' . (12) 
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Black solitary waves correspond to < u'^ < u 



2 ^ „,2 

oo' 



T = aV{u'J, /o = 0, So = a U{u^)-U{u)-ui,V{u'J , (13) 



and obey 

Finally, grey solitary waves correspond to < minjUQ, u"^} <u'^< maxJMQ, m^}, 

g {[[/(^o) - U{u^)] ul + (^L - uDViuDuU 
,2 _ 2a [[/(mqo) - Ujup) - {ul^ - ul)V{ul^)] m^m^ 

and are governed by 
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When Moo = or Wo = 0, one recovers from Eq. (TT6l) the bright and dark 
Eqs. (1121) and (TT4l) . respectively. Moreover, depending on the given forms for 
V{u?) and f/(-u), together with the particular values for u^o and Mq, the non- 
linearity coefficient a can be positive or negative. 



In the next sections we carry out the analysis for bright, dark and grey sohtary 
wave famihes when the nonhnearity is of the saturable form 




^K)=^^^ . (17) 



and therefore U{u) = — (1 + m^)^/(1 + u^). Notice that the cubic-like nonhn- 
earity (e.g. in Kerr media) corresponds for Eq. (1171) to the limit when m^ -C 1; 



that is, for V{u^) = (1 + m^)2(1 - 2u^) and U{u) = (1 + m^ 



,2^ _ f^ I „,2 \2f^ _ r,..2\ ^^A TTln.\ — l^ I „^2 \2/^2 _ ^4\ 



3 Qualitative Analysis 



Before deriving explicit localized solutions to Eqs. (IT^ . flT^ . and fITBl) when 
V{y?^ is given by flTTl) . we will carry out a qualitative analysis of Eq. (jH]) with 
the latter form for V"('U^): 

Let us first examine the case /o = 0. One easily verifies that Eq. ( ITSl) has three 
possible equilibrium points 



« = 0, Wi = ±W(l + n^)y^-l, (19) 



if (1 -|- M^) > A/r/o". When a < and F < 0, m = is a saddle point and 
u± are centers [see Fig. [U^a)], whereas for cr > and F > 0, m = is a center 
and u± are saddle points [see Fig. [l](b)]. In Fig. [I](a), the closed orbits of the 
phase-plane portrait {u,du/d^) correspond to periodic solutions of Eq. flTSj) . 
It is apparent that between the external and internal closed orbits there exists 
a homoclinic orbit (black curve), which, will be identified as a bright soliton 
solution in the next section. With respect to Fig. [Il^b), we have a heteroclinic 
orbit (black curve) in the phase-plane portrait, and, as we shall see in section 
5, it will correspond to a black soliton solution. Inside this orbit, there exist 
closed orbits which, again, represent periodic solutions of Eq. ( [T8l) . Outside 
the orbit, the solutions are no longer periodic. 

When /o 7^ 0, Eq. (ITS]) has a repulsive singularity and possess four nonzero 
equilibrium points Uj, j = 1,2,3,4; they correspond to two centers Ui = —U2 
(we omit their lengthy expressions), and two saddle points M3 = — M4 = Uoo- 
Two homoclinic orbits appear [see black curves in Fig. [U^c)], which will be 
interpreted in section 5 as grey solitary waves. Similarly with the cases in 
which /o = 0, the inner closed orbits Fig. [U^c) are periodic solutions. 
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Fig. 1. Phase-plane portrait of the real solutions to Eq. (jlSp for (a) a = 2T < and 
Uoo = 0; (b) (7 = r > and Uoo = 0.9; (a) and (b) both with /q = 0. (c) Iq / 0, 
a = 1.2r > and Uoo = 1-2. 

4 Bright Solitary Waves 



4-1 Explicit expressions for bright solitons 



Fundamental bright solitary waves fulfill Uqo = Uq = 0, and, without loss of 
generality, we may assume that uq > 0. Using Eqs. ( TTT1) . ( IT2l) . and (TT7I) . the 
constant F = cr/(l -|- Mq), and the resulting energy equation is 




2a (ul 



u^)u^ 



{l + ul){l+u^) 



(20) 



Now, since u{^) is a real function, it is thus necessary that a < 0. Equation fl20l) 
can be integrated exactly and it yields an implicit relation for the envelope 
distribution u{S,) of bright solitary waves 
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\ 2|a| 
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(21) 
Notice that, in this case, the obtained solution corresponds to the homoclinic 
orbit or loop of Fig. [T](a). Bright solitary profiles given by Eq. (!2T!) are shown 
in Fig. [2]^a). In the limit m < Mq -C 1 we recover from Eq. (1211) the well-known 



profile m(^) = ±Mosech['Uow2|o"| (^ — ^o)] of the one-dimensional bright soliton 
in a nonlinear cubic medium. Moreover, the dependence of the full width at 
half maximum (FWHM), A^, as a function of the peak amplitude Mq can be 
easily obtained from Eq. fl2T|) by setting u'^{A^/2) = Uq/2, and reads 



A^uo) 
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Figure EJ^b) displays the form of Eq. (l22l) . It is also known as the exis- 
tence curve for bright solitary waves because it provides a relation between 
the FWHM, the peak amplitude of the beam and the nonlinearity strength 
(through the parameter a) for which bright solitary waves can exist. Such a 




Fig. 2. (a) Transverse profiles of bright solitary waves for Uq = 3 and Uq = 0.5, 
as obtained from Eq. (|2ip . (b) Existence curve of bright solitary waves, where the 
FWHM A,^ is plotted versus the peak amplitude uq. 



dependence is often measured in experiments [32]. In our present case there 
are several interesting features. A property of A^(mo) is the so-called amplitude 
bistability. This means that for a given FWHM there are two different values 
of Uq for which bright solitary waves are to be found. This is depicted in Fig. 
I2]^a), where the two exhibited bright solitons are characterized by the same 
FWHM. It is evident that with increasing peak intensity, the bright solitary 
waves have a wider width. This is because the nonlinearity becomes weaker 
for large uq, and, therefore, in order to arrest the broadening (due to diffrac- 
tion), A^(mo) has to be larger. All these features are inherent of bright solitons 
in many saturable nonlinearities and constitute an ingredient for stability of 
the solutions [22]. The minimum of the existence curve can be evaluated by 
setting dAC,/duo = in Eq. fl22l) . and is obtained for uq — 1.059 (which gives 
3.338). These values are in exact agreement with those calculated 



numerically in Ref. [31] . 



4-2 Stability 



So far, we have been able to develop a full analytical approach allowing us to 
determine the exact profiles and the existence curve of bright solitary waves. 
However, an important aspect still remains open: Are the found bright solitary 
waves linearly stable with respect to small perturbations ? Using a well-known 
stability criterion for fundamental bright solitary waves [11132], we answer to 
this question. All we need to do is to calculate the dependence of the power Pg 
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Fig. 3. Dependence of the power Pg on the propagation constant, where 7 = r/|o"| 



of the bright soUtary waves on the propagation constant Fq. More explicitly, 
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where we have employed the symmetry property of the squared amplitude 
u^(— ^) = u^(0 and Eq. fl20l) . After integration, the result is 



Ps(%) 



1 + M^ 
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Using F = o"/(l + Mq), the power Pg can be cast in the form 
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(25) 



where 7 = F/|cr|. Fundamental bright solitary waves are stable ii dPs/dTo > 0, 
and unstable otherwise. Figure [3] shows that Pg is an increasing function of the 
propagation constant Fq for v'^/2 — \a\ < Fq < f^/2. Accordingly, the found 
bright solitary waves are all of them stable. 



5 Dark Solitary Waves 



5.1 Black Solitary Waves 



Fundamental black solitary waves fulfill mq = (mqo 7^ and u'q ^ 0). Using 
Eqs. (fT3l) . (IT^ . and (fTTl) we find that V = a and the energy equation is given 
by 
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2^2 



1 + m2 

Again, since n(^) is a real function, it is necessary now that cr > 0. 



(26) 



Equation (!26|) can be integrated exactly and it yields an implicit relation for 
the envelope distribution of black solitary waves 
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(27) 



The obtained solution corresponds to the heteroclinic orbit of Fig. [T](b). Black 
solitary profiles given by Eq. (l27|l are shown in Fig. |4](a). In the limit 1^1 < 
■Uoo ^ 1 we recover from Eq. (1271) the expression u{^) = ±'Uootanh['Uoo\/2a (^ — 
^o)] corresponding to the one- dimensional black soliton of the cubic NSE. 
Moreover, the dependence of the FWHM, A^, as a function of the amplitude at 
infinity Mqo can be easily obtained from Eq. fl27|) by setting m^(A,^/2) = m^/2, 
and reads as 



Ae(Moo) 
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As depicted in Fig. l4l(b), the amplitude bistability property of bright solitary 
waves is absent here. For increasing values of Uqo, the FWHM tends to the 



nonzero constant ln(2)/-i/2|o"|. 

Regarding the stability of the found black solitary waves, they are all stable 
since the saturable nonlinearity flT7|) belongs to the class of local nonlinearities 
for which the resulting solutions are stable for any velocity v [43]. 



5.2 Grey Solitary Waves 



Fundamental grey solitary waves fulfill uq ^ 0, Mqo 7^ 0, and u'q = 0. From 
Eqs. ( ITSal) . (ITSbj) . (lT6l), and i^ we find that F = a/(l + ul) and /g = 
(2crMoM^)/(l + Mq). This last expression imposes that cr > 0. The energy 
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Fig. 4. (a) Transverse profiles of black solitary waves for n^ = 2 and u^ = 1/2 
as obtained from Eq. ([27|) . (b) Existence curve of black solitary waves, where the 
FWHM A^ is plotted versus the amplitude at infinity Uoo- 



equation is given now by 




2a{u'-ul){ul-u') 

(1+m2)(1+m2)m2 



(29) 



Notice that the structure of Eq. fl29l) precludes the existence of solutions sat- 



isfying u^ < u < 



Un 



This means that there cannot be hump-like waves 



on a constant nonzero background; and so we look for dip-like waves (i.e. 



Wq < w < u^) in a constant nonzero background. Upon integration of Eq. (j29 
we derive the following implicit relation for the envelope distribution m(^) 
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(30) 



The obtained solution corresponds to the homoclinic orbits of Fig. [It^c). Fig- 
ure El^a) illustrates the profiles of two grey solitary waves having the same 
FWHM. The existence properties of grey-type solutions are described by 
an existence surface, which can be easily found from Eq. fl30l) by setting 



«2(A0 



Uif 



+ u^)/2. Figure Mjo) represents four slices of this surface for 



various ratios of Uq/uoo when varying u^. The amplitude bistability feature 
depends now on this ratio and disappears as Uq/u^o ~^ 0. Moreover, in contrast 
with bright and black solutions, the phase structure of grey solitary waves is 
nontrivial. The term j^ dT/v?{T) in Eq. ([6]) can be calculated exactly by re- 
sorting to Eq. (l29l) together with a method similar to the one employed to 
evaluate the power P^ of bright solitary waves [see Eq. (123|) ]. The resulting 
analytical form of the phase is 
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Fig. 5. (a) Transverse profiles of grey solitary waves for u^ = 4 and Uq = 1 (upper 
curve), and u'^ = 2.6 and Uq = 0.64 (lower curve), as obtained from Eq. ([30]) . (b) 
Existence curves for grey solitary waves, where the FWHM A^ is plotted versus the 
amplitude at infinity Uqo for uq/uoo = 0.1,0.4,0.7, and 0.9 (lower to upper curves). 
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where u = u{^) is given implicitly by Eq. fl30|) . In the limit Uq < ^^ < ^L ^ I5 
one retrieves from Eq. fl30|) the profile for the one- dimensional grey soliton in 
a nonlinear cubic medium, which reads as 



«(0 



±w«^ - («^ 



Mq) sech^ 



\/2R^4^^(e-eo) 



(32) 



and the corresponding phase dependence follows straight away from Eq. (131]) . 

With respect to the the stability of the found grey solitary waves, it is an 
open problem to determine the regions defined by Uq and Mqo for which these 
solutions are stable. This issue will be addressed in a subsequent work. 



6 Conclusions 



In this paper we have constructed explicit analytical expressions for solitary 
wave solutions of the nonlinear Schrodinger equation with a saturable non- 
linearity of the form oc (1 + jAp)"^ arising, for instance, in the propagation 
of nonlinear optical beams in centrosymmetric photorefractive materials. The 
bright and black soliton solutions had only been found numerically previously 
in the context of optical applications of the model equation. In the case of 
bright solitons solutions we have also performed a study of their stability. The 
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analysis of the stability of grey solitons is not so simple and could be the goal 
of future research in this area. 

Our work can be extended in several directions. The first one, by consider- 
ing more complicated saturable nonlinearities such as those arising in media 
with electromagnetically induced transparency [H]. Secondly, by studying sat- 
urable vector media with several fields involved. Finally, another interesting 
extension could be the analysis of the model with space-dependent parameters, 
resorting to the technique of Lie symmetries 
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